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• 1— I , Abstract 

;h ' A determination of the hadronic fragmentation functions of the boson is 

^ presented from a study of the inclusive hadron production wit the DELPHI detector 

at LEP. 

These fragmentation functions were compared with the ones at lower energies, 
thus covering data in a large kinematic range: 196 < < 8312 GeV^ and 
x{= Ph/ Ebeam) > 0.08. A large scaling violation was observed, which was used to 
extract the strong coupling constant in second order QCD: 

as{Mz) = 0.118 ±0.005. 

(5) 

The corresponding QCD scale for five quark fiavours is: ^-jj^ = 230 it 60 MeV. 



(Contribution to the XVI International Symposium on 
Lepton-Photon Interactions, 
CorneU, 10-15 August, 1993) 
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1 Introduction 



Hadron production in e^e~ annihilation originates from the production of quark-antiquark 
pairs, which can radiate gluons, the quanta of the field theory of the strong interactions, 
Quantum ChromoDynamics (QCD). Gluon radiation depends logarithmically on the cen- 
tre of mass energy due to the increasing phase space with increasing energy and the energy 
dependence of the running coupling constant of QCD. These effects lead to variations of 
the momentum spectra of the produced hadrons as a function of the centre of mass energy, 
even if the momenta are scaled to that energy. These scaling violations can be used to 
determine the strong coupling constant a^- 

For example, the scaling violation in deep inelastic lepton-nucleon scattering leads to 
as{Mz) = 0.112 ± 0.005 ll|, |, H. This is somewhat lower than, but not in disagreement 
with, as measurements at the Z° mass from shape variables, jet rates and total cross 
sections as measured at the electron-positron storage ring LEP|^, |^. 

Until now has not been determined from scaling violations in e"'"e~ annihilation, 
since this requires precise data at two very different energies to observe a significant 
scaling violation. Data collected at the PEP and PETRA storage rings were only precise 
at energies around y/s = 30 GeV@, ^ |[ and the scaling violation was only observed 
qualitatively 0. 

In this paper we present data of the inclusive hadron spectra, as measured with the 



DELPHI detector []I0[ at LEP and present the first determination from the scaling 
violations in the fragmentation function by combining the data at LEP with data from 
the PEP, PETRA and TRISTAN storage rings. The squared four-momentum transfer 
from the incoming leptons to outgoing hadrons studied here is two orders of magnitude 
larger than the ones studied in deep inelastic scattering, thus avoiding regions where 
non-perturbative effects noticeably influence the results. On the other hand, in e~^e~ 
annihilation one has to combine data from different experiments at different energies and 
study the effect of the varying quark flavour composition due to the Z^-resonance. As 



will be shown, these are not dominant uncertainties (see also Ref. 11 



2 Determination of the Fragmentation Function 

The inclusive production of charged hadrons in the reaction e+e" ^ h + X can be 
described by two kinematic variables, and x, where is defined as the square of the 
four-momentum transferred from the leptons to the hadrons and x is the fraction of the 
beam energy Ebeam carried by the hadron h. In e~^e~ annihilation equals s, the total 
centre of mass energy squared. 

The fragmentation function D{x, Q"^) is directly related to the scaled hadron momen- 
tum distribution: 

D{x, Q') ^ W,{Q')Di{x, Q') ^ -^(e+e- ^h + X), (1) 
~i (Tt ax 

where at is the total cross section and D{x, Q"^) is the sum of fragmentation functions 
Di{x,Q'^) over all five flavours, each having a weight Wi{Q'^) given by the electroweak 
theory. 

The DELPHI data were collected during 1991 at energies near the Z° peak. Multi- 
hadronic events were selected according to the criteria given in Ref. ||12||. The selection 
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required that there were at least 5 charged particles with momenta above 0.2 GeV/c and a 
track length in the detector of at least 50 cm, that the total energy of the charged particles 
exceeded 15 GeV (pion mass assumed), that the forward and backward hemisphere with 
respect to the beam axis each contained a total energy of charged particles larger than 
3 GeV, and that the polar angle of the sphericity axis was between 40" and 140°. In 
addition the momentum imbalance was restricted by requiring that the absolute sum of 
the three-momenta of charged particles was less than 20 GeV/c. 

After these selection criteria, 186774 events at a mean centre of mass energy of 91.2 
GeV were kept. The background due to beam-gas scattering and 77- interactions was less 
than 0.1% and r+r^ events contributed 0.2% to the selected sample. 

The scaled inclusive momentum spectrum was obtained by correcting the x distribu- 
tion of all charged particles for initial state radiation, particle decays, detector effects, 
and selection cuts. In principle x is the fraction of the beam energy Ebeam carried by 
the hadron h, i.e. x = Eh/Ebeam, but instead of Eh the momentum ph was used. This 
was experimentally better measured and provided the same scaling violation information. 
The corrections were obtained from a detailed Monte Carlo simulation of the detector 
after generating the hadronic decays with the Monte Carlo program from the Lund 



group |]T3|. Higher order initial state bremsstrahlung radiation was implemented by using 



the DYMU3 program]!^. The corrected data were obtained by multiplying the data in 



each bin of a histogram by a correction factor defined as 



(i) _ ^"gen /px 



gen 
Sim 

where A^^eli contents of the histogram bin i at the generator level witout initial 

state radiation, and Ng2n after initial state radiation and detector simulation. All primary 
particles with a lifetime larger than 3 ■ 10~^° s were assumed to be stable at the generator 
level, i.e. they were included in Nj^^}^, and all those with a shorter mean life (including 
Kg^s and A's) were allowed to decay as part of the simulation process. Thus the corrected 
distributions include the contributions of these short-lived particles, as was the practice 
in experiments at lower energies. Each histogram was normalized to the total number of 
events. The corrected distribution is simply: 

^^corr — ^ Data- K'^) 

The corrections varied smoothly and were less than 35% for x below 0.8 (see Fig. ^a). 
Only this range was used for the comparison with QCD in the next section. The correction 
factors deviated from one due to acceptance losses and momentum measurement errors. 
For larger x values the momentum measurement errors dominated and, together with 
the steeply falling spectrum, caused a smearing towards higher momenta, resulting in a 
correction factor of 0.4 near x=l. The corrected spectrum and thus the fragmentation 
function was obtained from the total number of events A^^^ and the corrected event numbers 
^corr Qd-ch. X value: 

D{x, Q') ^ ^ ^h + X), (4) 

i\t d,x at ax 

The corrected spectrum is displayed in Fig. ^b and tabulated in Table |l| together with 
the statistical and systematic uncertainties. 
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The systematic uncertainties were determined by varying the selection criteria and 
by using different Monte Carlo simulations. The largest uncertainty in the correction 
factor is connected with the uncertainty in the charged multiplicity: the integral of the x 
distribution is equal to the total charged multiplicity, since each event has Nch entries, so a 
systematic error in the multiplicity after detector simulation translates into an uncertainty 
on the normalisation of the x distribution. Varying the cuts, especially changing the 
minimum number of charged particles from 5 to 6 and varying the cut on the sphericity 
axis between 30° and 45°, changed the correction factors by less than 10% of their deviation 
from 1, i.e. typically 3% in the intermediate range, but up to 10% for x above 0.8 and x 
below 0.04. 

In the LUND Monte Carlo program several generators can be used. Partons can be 
generated either with the "Parton Shower" algorithm or by using the exact O(as^) QCD 
matrix element. The difference in correction factor between these two options was less 
than 1% in the intermediate x range. More details can be found in Ref. [1^ . 



The relative systematic uncertainty from the sources mentioned above was parametrised 
in the following way: 

5^1 = maa;(0.03, 0.1 | 1 - c^'^ |) (5) 

This procedure gives a relative error of at least 3% for the intermediate x range and in- 
creases the error near the endpoints. It should be noted that the systematic uncertainties 
are correlated between the bins, since a change in the selection criteria moved the correc- 
tion factors for each bin all in the same direction. These correlations will be taken into 
account in the determination of the strong coupling constant. 

A significant scaling violation is observed between the DELPHI spectrum at 91 GeV 
and the data from TASSO0 at a lower centre of mass energy (see Fig. [^). Note that the 
errors are smaller than the symbols for most of the data points on this logarithmic plot. 
In order to show the scaling violation and the errors more clearly, the ratio of the curves 
in Fig. |l]b is shown in Fig. ^. As can be seen, the scaling violation varies from + 40 % to 
-30% for X varying between 0.01 and 0.7. The deviation of the data from the horizontal 
line in each x bin is, to first order, proportional to the strong coupling constant, so each 
data point yields an independent measurement of as- All values should be consistent, 
which is a strong constraint and simultaneously a cross check. 

In Fig. |]the DELPHI data are compared with other experiments at lower energies|^ 
0, I], |T6| for several x intervals. Clearly, at small x values the fragmentation function 
increases about 30%, while at high x values it decreases by roughly the same amount. 
This is exactly what is expected: the higher the energy, the more phase space becomes 
available for gluon radiation. Since the primary quarks lose more energy due to radiation, 
this depopulates the high x region. The radiated gluons tend to populate the small x 
region, increasing the spectrum there. 

The curves are the QCD fits for high and large x, as will be described in the next 
chapter. Clearly all x ranges agree well with the QCD fits, even if they are extrapolated 
to small X and small Q^. The fact that all regions can be described by a single value of 
the QCD scale A^g^ provides the cross check mentioned above. 



3 Comparison with QCD 
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3.1 Theoretical Framework 



The scaling violations in the fragmentation function, defined by Eq. 0, are described by 
the coupled integro-differential evolution equations[0], which can be written as: 



dQ^ \ Dg{x,Q^) / 27r 



(6) 



V 



The splitting functions Pij{z) in the 2x2 matrix are the probabilities of finding parton 
i with momentum fraction z from its parent parton j where i,j = G refers to a gluon 
and = g to a quark. Note that a gluon can split into a quark-ant iquark pair of any 
flavour, hence the summation and the factor two in front of Pqc- As mentioned before, 
D represents the sum over the weighted contribution of each flavour (see Eq. |I|), each 
having its own fragmentation function Di. 

The splitting functions can be obtained by integrating the exact QCD matrix element. 
In order to obtain the probability of finding a hadron with momentum fraction x, one 
has to integrate Pij{z) convoluted with the probability Di{x/ z^Q"^) that the parton with 
energy fraction z fragments into a hadron with momentum fraction x. The symbol ® 
denotes a convolution integral: 

P,,{z) ® D{x, Q') = / -P,,iz) ■ D{-, Q'). (7) 
J z z 

X 

Note that x/ z is the fractional hadron energy, i.e. xjz = Ph/Vvarton, since x = Ph/Ebeam 
and z = Pparton/Ebeam- Obviously, z has to be larger than x, hence the lower bound in the 
integral. 

The evolution equations describe the dependence of the fragmentation function. 
Their solutions have not yet been found in an analytical form. Numerical solutions, which 
account for second order corrections to the splitting functions or to the anomalous dimen- 
sions have been developed in Ref. Alternatively, one can integrate the exact second 



order QCD matrix element directly, which has some advantages, as will be discussed later. 

In principle the fragmentation of quarks involves an infinite number of soft and 
coUinear gluons. Hence, a cut-off on the isolation of the gluons is used in order to de- 
cide whether a gluon should be part of the quark fragmentation or if it should fragment 
independently. In the latter case it contributes to Dq instead of Dg. 

As a cut-off, the minimum invariant mass between quarks and gluons was required to 
be above 9.1 GeV/c^. This cut presents an arbitrary definition of quarks and gluons, but 
it has to be made in any analysis of the scaling violations. It was varied in order to study 
its effect on the determination of A^^, as will be discussed in the section on systematc 
errors. 

Such an invariant mass cut selects a certain part of phase space, which varies with 
energy, as shown in Fig. ^. Here the energy fractions Xk = Eq/Eheam of both quarks in 
second order QCD were plotted against each other at centre of mass energies of 35 and 
91 GeV. The 2-jet events are located at Xi = X2 = 1 and the 3- and 4-jet events more 
towards the centre. One clearly observes the strong increase in phase space for the events 
away from the 2-jet region. The invariant mass cut 

y = M^/s = 1 - Xfc, 
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where Mij are the invariant masses between any pair of partons and Xk are the fractional 
quark energies, ehminates the soft and colhnear gluons in the regions 1 — Xk < Mf-/ s = 
0.0676(0.01) for the centre of masss energies of 35 (91) GeV. The difference in phase 
space between these energies increases the qqG cross section with a given 'hardness' of 
the gluon, i.e. with a given invariant mass cut, by a factor four: the 3-jet rate varies from 
20% to 80%, as shown (for a constant value of as of 0.121) by the solid line in Fig. | . 

In addition to the dependence of the phase space, one has to consider the 
dependence of as, which has the opposite effect: it decreases the qqG rate with increasing 
energy. This decrease, from the running of the coupling constant, can be observed if the 
phase space for the qqG final state is defined as a constant fraction of the total phase 
space, for example by a constant y-cut instead of a constant invariant mass cut. In this 
case the only dependence comes from the running of the coupling constant, which 
decreases the qqG cross section by about 20% if the centre of mass energy is increased 
from 35 to 91 GeV[Q. This decrease of the 3-jet rate in a constant fraction of phase space 
is shown (for a fixed value of of 270 MeV and renormalisation scale = s) as the 

^ MS ^ ' 

dashed line in Fig. |g. Note that the scaling violation from the running of as is a small 
effect compared with the scaling violation from the change in phase space, as is apparent 
from Fig. |^. 

The large phase space dependence can be absorbed in the fragmentation function, 
which then depends on both x and Q^. The redefined cross section has a well determined 
perturbative expansion in ^^(Q^). This would not be the case if the large phase space 
corrections, proportional to aglnQ^, were considered to be QCD corrections. 

The energy dependence of as can be expressed in terms of the energy independent 

(5) 

QCD scale A^^; here the upper index indicates the number of flavours nf = 5 and the 
lower index the renormalisation scheme (following the convention of Ref . |]19| ) : 

4:71 



as{ii 

with 



^ Pi InL 



PS L 



(8) 



(3o = 11 

A = 2(51 -fn;) 
The energy scale /x^ of as can be related to Q'^ = s by 

/^^ = fs, 

where / is the renormalisation scale factor. The choice of / is free and QCD predictions 
would not depend on it if all higher orders were known. In practice, calculations have been 
performed only up to a fixed order and varying / in a wide range indicates the uncertainty 
due to the higher orders, as will be discussed in the section on systematic errors. Note 
that different choices of / change the value of ag- In order to keep the physical observables 
constant, the coefficients of the higher order terms in the expansion of the observable 
have to be changed correspondingly^. 

The extraction of as from the observed scaling violations is straightforward. First, the 
X dependence of the fragmentation function, which cannot be calculated perturbatively, 
must be parametrised from data at a reference energy. Starting from this parametrisation 
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at the reference energy, the evolution to higher energies is predicted by QCD and compared 
with the observed fragmentation function at these energies. In the following sections 
the parametrisation of the x dependence and the dependence of the fragmentation 
functions will be discussed. 



3.2 Parametrisation of the Fragmentation Function 

The fragmentation functions have been studied in great detail in e~^e~ annihilation. Even 
such details as the "string effect", predicted in QCD by the interference effects of multiple 
gluon emission, have been observed] 2^] and can be well described by the string fragmen- 



tation model[0. Although any parametrisation of the x dependence at a given would 
suffice, we have chosen the string model for the following reasons: 

• The quark and gluon fragmentation functions are described by the same string with 
the same parameters, thus reducing the number of free parameters. 

• In this model soft gluons are automatically "absorbed" in the string, i.e. they only 
produce some transverse momentum to the string, but do not lead to independent 
jets. In independent fragmentation models the fragmentation of soft gluons is prob- 
lematic because of phase space restrictions for hadron production. 

• Quark mass effects are taken into account in the string model. 

Hadrons inside a jet are characterised by the limited transverse momenta with respect 
to the jet axis independent of the jet energy and the longitudinal momentum spectra. 
These momentum components can be parametrised by two energy independent functions, 
a longitudinal and a transverse fragmentation function. Italics have been used here in 
order to distinguish these parametrisations at a reference energy from the fragmentation 
function D{x, Q"^). 

The transverse momentum spectrum was parametrised by a Gaussian with a variance 
of (500 MeV/c)^|21|. The longitudinal momentum spectra of light and heavy quarks are 



parametrised differently since the latter have much harder spectra because of their larger 



mass. The Lund symmetric fragmentation function |22| was used for the light quarks: 



f(y) = y — j ' (9) 

where mj_ = \J m? + p]_ is the transverse mass of the hadron, y = {E + pi)^/ {E + pi)p 
determines the fraction of the primordial parton energy taken by the hadron h, with p, 
indicating the parton with energy E and longitudinal momentum pf, a and b are two free 
parameters. The longitudinal spectrum depends mainly on a — 6 which scales like N, 
where N is the total multiplicity, so effectively there is only one free parameter. For the 
heavy quarks (6 and c) the Peterson fragmentation function||2^ was used: 



m = - 

y 



i-i- 



y 1 



(10) 



Here the free parameter, e^, is expected to vary as 1/m^, so ec/e?, = 9.4 was used. Hence 
there are only two free parameters to tune the momentum spectrum (a and e?,). The pa- 
rameters determining the fraction of strange quarks picked up from the vacuum, the ratio 
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of vector to pseudoscalar mesons, the fraction of baryons, as well as the decay parameters 
were all left at their default values, since a good description of the x dependence was 
possible with these. 



3.3 Determination of the Strong Coupling Constant 

The dependence of the fragmentation function can be derived either from the evolu- 
tion equations or from a direct numerical integration of the exact QCD matrix element. 
Since the splitting functions have been derived from the integration of the matrix element, 
both methods are, in principle, equivalent. However, higher order differences might oc- 
cur. As will be discussed in the section on systematic errors, these differences are small. 
Therefore, the dependence has been determined from the integration of the exact 



second order QCD matrix element, using the formulae from reference P3, which have 



been implemented in the Lund string model|T^. This method has the advantage that the 



convolution of the splitting and fragmentation functions is done in a consistent manner, 
i.e. the cuts to separate the nonperturbative region from the perturbative one are the 
same for the splitting and fragmentation functions. These cuts will be discussed in more 
detail in the section on systematic errors. Furthermore, the weights Wi{Q'^) in Eq. |l|from 
the electroweak theory have been incorporated in this model. 

The strong coupling constant was extracted in the following way. A simultaneous fit of 
the QCD scale and the fragmentation parameters a and e^, was made by minimising: 



(f^Di^-T(^)A (1 - /,) 



(fjCTexp) 



(11) 



where fj is the normalisation factor for experiment j with data D^^^ in a given x bin with an 
experimental error a^^^p for that bin and an overall normalisation error cr„. The fit function 
was the x parametrisation from the string model convoluted with the dependence 
from the integration of the exact QCD matrix element and the dependence of as- In 
order to prevent a bias from b-quark threshold corrections, only data at or above i?cm=29 
GeV were used . Furthermore, data at high and low x values have not been used, since the 
experimental correction factors are large in these regions. Fitting the data from Delphi 
simultaneously with all other available data|^, [7|, |^, ^ in the range 0.18 < x < 0.8 



and < Q2 < 9X 22 Q^y2 yielded the results given in Table 

The fit was repeated for two values of the renormalisation scale. For / = /x^/s = 
0.01(1.0) the result was: A^^ = 193t?? {^QdtH) MeV. The fit resuhs were obtained for 
a value of 6 = 0.283 in the Lund symmetric fragmentation function^. A good agreement 
was observed for all x values with the same fragmentation parameters at both 35 and 
91 GeV, so the difference between the energies depended on A-^ only. The results for 
f = 0.01 are shown as the solid lines in Fig. |I|b; the x of the fit for / = 1.0 was equally 
good. 



""^ Although a and b are strongly correlated, one could not leave b at an arbitrary value and just fit a 
or vice- versa. A good parametrisation was obtained if b was chosen in a range around 0.3. 
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3.4 Systematic Uncertainties 

The results in the previous section include both systematic and statistical uncertainties, 
as well as the uncertainties from the correlation between the fragmentation parameters 
and A^^. In addition there are theoretical uncertainties from the unknown higher order 

MS 

corrections, which are usually estimated by varying the renormalisation scale. To get the 
complete error estimate, the following have been investigated: 

• Experimental uncertainties. 

In the definition of x^? E<1- ^ixp represents the total error for that data point, 
obtained by adding in quadrature the statistical and point-to-point systematic un- 
certainty, but excluding the overall normalisation error, o"„. However, the separation 
between point-to-point systematic uncertainty and normalisation uncertainty is not 
straight forward and usually not given in the literature. Furthermore, the published 
systematic uncertainties are not always comparable in the possible sources which 
have been included. Therefore the systematic uncertainties were varied considerably 
in order to check their influence on the fitted value of A^^. The following procedure 
was adopted: the total error, atot, was split into a point-to-point error a^l^p and a 
normalisation error (T„: 

= ^S/ + (12) 

(Tn was varied from 1% to 3% and subtracted from the total error quoted by the 
experiments (using Eq. 0). If the remaining point-to-point error fell below a certain 
minimum value, it was adjusted to this minimum value in order to ensure that the 
point-to-point error squared would not become negative or too small for experiments 
in which all possible systematic effects had not been included in the error estimate. 
This minimum value o"^™ was varied between 1% and 3%. Of course, the of the 

(5) 

fit was changed if the errors were changed, but fortunately A^g^ varied by only ±3 
MeV if (Tn and 0"™^" were changed in the ranges given above. These small changes 

in A^^ indicate that the scaling violations are determined only by the shape of the 
distributions, not the absolute normalisation. The values have been summarized 
in the Tables ^ using a„=2% and al^'J^ = 3%. The total is 73 for 71 data 
points in the fit region using a renormalisation scale factor of 0.01. For comparison, 
the values outside the fit range are shown too. The correlations between the 
parameters depended on the assumed errors, but were never larger than 40% for 
any pair of parameters. The fit normalisation factors were consistent with one for 
all experiments, as indicated in the captions of Tables |T| and ^-|^. 

• Differences between experiments. 

In contrast to the deep inelastic lepton scattering experiments, which measure the 
dependence in a single experiment, we had to combine data from different ac- 
celerators, so systematic effects from differences between experiments had to be 
considered. They were checked by comparing the results of different combinations 
of experiments. The maximum difference in A^^ from the various combinations 
of the 6 experiments was less than 30 MeV, which is not much larger than the 
statistical uncertainty (see Fig. P). The systematic uncertainty from this source 
was conservatively estimated to be half the maximum difference, i.e. 15 MeV, thus 
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assuming that the whole difference was systematic and not due to statistical fluc- 
tuations. The reason for this surprisingly small spread is simple: all experiments 
used large An solenoidal detectors in which the momentum spectrum, especially the 
shape, was clean and easily measured. As mentioned before, it is the change in 
the shape of the x spectra which determines the scaling violation, not the absolute 
normalisation. 

x-dependence. 



For low X values the contributions from multiple soft gluon emission start to domi- 
nate. In this region the of the string model parametrisation becomes somewhat 
worse (see Tables |l] and To estimate the uncertainty from the small x range, 
we fit between x^m and Xmax and varied Xmin between 0.08 and 0.4. Since the ex- 
perimental correction factor for high momentum particles becomes large for x > 0.8, 
Xmax was kept at 0.8. For Xmin =0.08, A^g^ for / = 0.01 increased from 190 to 210 
MeV, but for the fit range considered {xmin > 0.18), no variation in as was seen 
(see Fig. Nevertheless, the uncertainty for the selected x range was conser- 
vatively estimated to be 10 MeV, which is half the difference between the values 
obtained for Xmin = 0.08 and Xmin =0.18. As mentioned in the introduction, each 
x- value provides an independent determination of a^. The fact that is practically 
independent of the selected x-range indicates that all x- values are consistent. 

Heavy Quark Fractions. 

The fragmentation effects largely cancel in the difference between the spectra at 
different energies. However, since the primary quark composition changes with 
energy, the influence of the difference in fragmentation between light and heavy 
quarks should be considered. Although the primary mesons from heavy quarks have 
the hardest momentum spectra, the spectra after decays are not much different from 
the ones for the light quarks and actually somewhat softer. Furthermore, it was not 
possible to mimic the characteristic change in shape from the QCD scaling violations 
by the difference in quark compositions, as shown by the dashed-dotted line in Fig. 



Fitting the x spectra at 35 and 91 GeV simultaneously was a good way to determine 
the fragmentation of both light and heavy quarks, since the different quark com- 
positions at the different energies, combined with the somewhat softer x spectrum 
of the heavy quarks, yielded only a moderate correlation between the fragmenta- 
tion parameters a and (see Table The fitted value of the latter parameter 
gave an average x of the B-hadrons of 0.69 ± 0.01 at the LEP energy, which is 
in good agreement with the value obtained from lepton spectra in semi-leptonic B 
decays p6|, p7| . Note that the determination of ef, from the inclusive hadron spectra 



included all decays and was therefore independent of the value determined from the 
lepton spectra. 

As an additional check that the different heavy quark fractions at 35 and 91 GeV do 
not mask the scaling violation from QCD, the scaling violation was calculated with 
a constant fraction of heavy quarks (~ 11% for b-quarks and ~ 44% for c-quarks, 
which are the values at ^/s=35 GeV). The amount of scaling violation is not changed 
significantly, as shown in Fig. the small difference was taken into account in the 
fit and the residual uncertainty in A^^ was estimated to be 10 MeV. 
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Independent versus String Fragmentation. 

In the fit, the string fragmentation model was used to parametrise the x depen- 
dence. As an alternative, the independent fragmentation option in the Monte Carlo 
program from the Lund group has been usedQ. In this case all quarks and gluons 
fragment independently. The whole analysis, including the parametrisation of the 
X dependence, was repeated with this model. The fit quality was similar and the 
value of Aj-j^ was not changed outside the experimental errors, again indicating that 
fragmentation uncertainties largely cancel in the difference between the spectra at 
different energies. Half the difference between the different fragmentation models 
(9 MeV) was taken conservatively as the error for fragmentation. 

The renormalisation scale uncertainty. 

(5) 

As mentioned before, the renormalisation scale is a free parameter; A.-^ would be 
independent of the choice of this scale if all higher order corrections were known. 
However, in a fixed order calculation a lower scale implies a larger value of as- For 
the 3-jet cross section the change in the Born cross section can be compensated 
by a different coefficient in the higher order correction. However, the 4-jet cross 
section is only known up to the Born term in second order QCD, so a lower scale 
for the argument of implies a higher 4-jet rate. The of the fit did not change 

(5) 

significantly by changing the scale, but varied from 190 to 270 MeV if the scale 
was changed from to Ecm, which corresponded to a change in from 0.115 

to 0.121 (see Fig. |^). Thus this error, originating from the unknown higher order 
corrections, has been found to be dominant, as in all other determinations]^]. The 
scale dependence was still relatively small, since we studied the difference between 
the spectra at different energies, so higher order contributions and fragmentation 
effects partially cancel. Similar observations hold for other "difference" variables, 
like the Asymmetry in the Energy-Energy Correlations (AEEC) or the difference in 
jet masses [||. 

An independent estimate of the higher order contributions can be obtained from the 
equations (^: in these equations the higher order terms are taken into account by 
exponentiating the leading logarithms proportional to (as/27r)" In" Q^, which appear 
as leading terms in a calculation to order n. The difference of these terms between 
Q = Qmin and Q = Qmax IS proportioual to (as/Tr)" ln" (Qmax/Qmm)- Since in our 
case the difference in Qmax and Qmin is only a factor three, the exponentiated form 
of the leading logarithms will be close to its second order expansion, so the higher 
order contributions are expected to be small. This can be checked explicitly by 
integrating these equations in n steps. Since at each step a gluon can be emitted, this 
corresponds to summing up all higher order terms proportional to (as/vr)" In" Q^. 
The change in scaling violation between 30 and 90 GeV was found to be less than 
5% if n was varied between 2 and 20, so this change is similar to the uncertainty 
from the scale dependence. Since this exercise was done only in first order, using 



the program from reference [28|, the larger range from the scale dependence was 



used as an estimate of the error from the unknown higher order corrections. 



^The following parameters were used: f=0.01,b=0.283, pt=500 MeV/c, stopping point of fragmen- 
tation 3.1 GeV/c. The gluon was splitted into light quark pairs according to the Altarelli-Parisi 
splitting function and each quark fragmented then independently. The fitted parameter values were: 



a = 0.88 ± 0.05, ef, = 0.007 ± 0.001, Atto = 211^36 MeV. 
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• Cut-off dependence. 

As mentioned before, the fragmentation of quarks involves a large number of soft 
and collinear gluons. Hence, a cut-off on the isolation of the gluons was used in order 
to decide whether a gluon should be part of the quark fragmentation or if it should 
fragment independently. Below this cut the quarks and gluons were considered to 
fragment into a single jet, i.e. they were considered to belong to the non-perturbative 
regime in the model and were recombined beforehand. 

As a cut-off, the minimum invariant mass between quarks and gluons was required 
to be above 9.1 GeV/c^. The scaling violations were not very sensitive to this cut, 
since they just required a different parametrisation of the nonperturbative part for 
a different cut. What mattered was a good parametrisation of the x dependence. 
The cut could not be decreased, since with this cut practically all phase space was 
already used at the highest energy, as shown in Fig. |[ Decreasing the cut further 
would cause the 4-jet cross section to become so large and positive, that the 3-jet 
cross section would become negative in some regions of phase space due to the large 
and negative virtual corrections in the second order QCD matrix element in that 
case. Increasing the minimum invariant mass squared by a factor two resulted in an 
mcrease of Ag^ of 60 MeV. Therefore an error of ±30 MeV was attributed, although 
part of this was presumably already absorbed in the scale error: increasing the cut- 

(5) 

off or increasing the renormalisation scale increased A^g^ in both expected 
for a decrease of the higher order contributions from multiple gluon radiation in 
both cases. This can be seen in Fig. H, which shows the phase space of 4-jet events 
at an energy of 91 GeV for the default parameters (/=0.01 and |/min=0.01) at the 
top and the changes if / is increased to 1 and ymin to 0.02. For these cases the 
total 4-jet rate decreases from 28% at the default value to 12 and 9%, respectively. 
Note that the change from the change in the ymin cut is not only at the border, 
as one might naively expect. This can be understood, if one realizes that the cut 
converts many 4-jet events into 3-jet events, which can be located anywhere in the 
allowed phase space. Consequently, the variations of i/min and / change the higher 
order corrections in a vary similar manner and by counting them as independent 
contributions to the total systematic error is a conservative attitude. 

The total errors were obtained by adding in quadrature the errors from the fit {lili MeV), 
from the x-dependence (10 MeV), from heavy quark fractions (10 MeV), from fragmen- 
tation (9 MeV), from the comparison between experiments (15 MeV), from the gluon 
cut-off dependence (30 MeV) and from the scale dependence (40 MeV). A summary of 
the systematic errors is given in Table ^ 

4 Summary 

A precise determination of the fragmentation function in the decay of the Z° boson has 
been presented. A comparison with the fragmentation functions at lower energies shows 
a strong scaling violation, which leads in second order to a QCD scale A^^ between 190 
and 270 MeV (see Table 0). Taking the average as the central value and using the total 
uncertainties as given in Table |g resulted in Ai^ = 230 ± 60 MeV, which corresponds to 

as{Mz) = 0.118 ±0.005. 
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These results in the time-hke region are in good agreement with the results on scaling 
violation from deep inelastic lepton nucleon scattering (space-like region; = 0.112 ± 
0.005 [H, 0, ^) and with other as determinations at LEP from jet rates and shape variables 
(a, = 0.120 ±0.007 gi). 

It should be noted that future e^e~ linear colliders can operate over a large energy- 
range, which allows a study of the scaling violations with a single detector. For example, 
a 500 GeV linear collider can measure a large scaling violation by running for a short time 
ath the 91 GeV Z° peak. From Fig. |^ it can be seen that the expected scaling violation 
will be even larger than the one observed between 35 and 91 GeV, thus allowing a precise 
determination of at high within a single experiment. 
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0.50 


0.561 
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0.396 
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0.376 


0.60 - 0.66 


0.66 


0.268 


0.007 


0.008 


0.260 


0.66 - 0.72 


1.97 


0.160 


0.006 


0.005 


0.167 


0.72 - 0.78 


3.53 


0.096 


0.004 
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0.78 - 0.84 


15.69 


0.045 


0.002 


0.002 


0.055 


0.84 - 0.90 


5.72 


0.023 


0.002 


0.002 


0.019 


0.90 - 1.00 


7.58 


0.0059 


0.0005 


0.0005 


0.0044 



Table 1: The inclusive hadron x spectrum as measured by DELPHI as well as the 
statistical and systematic errors. The prediction from the exact QCD Matrix Element 
calculation followed by string fragmentation is shown under the label QCD+SF and the 
of each bin is shown in the second column. The centre of mass energy is 91.2 GeV and 
the overall normalisation factor from the fit is 0.995 (not included in the data column). 
Only the data between the empty rows was used for the determination of ccg. 



17 





/ = 0.01 


/= 1.0 


MS 

a 


193+?? MeV 
0.85 ± 0.03 


269+11 MeV 

— 1^ 

0.96 ±0.03 




0.009 ±0.002 


0.008 ±0.002 


X^/datapoints 


1.02 


0.99 


corr. a — A-j4tt 


-0.07 


-0.06 


corr. eb - A)^ 


-0.36 


-0.31 


corr. eb — a 


-0.07 


-0.22 



Table 2: Results of the fit to all data (71 data points for > 29^ GeV^ and 0.18 < x < 
0.8) for two renormalisation scales (/ = n'^/s = 0.01 and 1.0, respectively). 
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3.75 
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2.12 
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0.023 
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0.70 - 0.80 


5.46 


0.162 


0.018 


0.119 




0.80 - 1.00 


0.10 


0.0300 


0.0120 


0.0259 



Table 3: As Table 1, but for the TASSO experiment [6]. The centre of mass energy is 35 
GeV and the overall normahsation factor from the fit is 0.984 (not included in the data 
column). The fourth column includes both the statistical and systematic errors. 
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0.008 
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Table 4: As Table 1, but for the CELLO experiment [7]. The centre of mass energy is 35 
GeV and the overall normalisation factor from the fit is 1.007 (not included in the data 
column) . 
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x-Bin 




Data 




QCD+SF 


0.00 - 0.05 


1.71 


115.4 


1.8 


115.7 


0.05 - 0.10 


3.26 


65.1 


1.2 


71.6 


0.10-0.15 


6.29 


31.6 


0.6 


30.5 
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6.33 


17.5 
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16.9 
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2.60 


10.40 


0.21 
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0.00 


6.29 


0.13 


6.58 
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0.07 


4.07 


0.09 


4.21 


0.35 - 0.40 


0.05 


2.76 


0.07 


2.90 


0.40 - 0.45 


2.01 


1.80 


0.06 


1.96 


0.45 - 0.50 


5.16 


1.18 


0.04 


1.32 


0.50 - 0.55 


0.37 


0.810 


0.039 


0.868 


0.55 - 0.60 


1.43 


0.515 


0.031 


0.574 


0.60 - 0.65 


0.23 


0.347 


0.023 


0.351 


0.65 - 0.70 


2.45 


0.227 


0.020 


0.269 


0.70 - 0.75 


0.43 


0.167 


0.020 


0.161 


0.75 - 0.80 


1.24 


0.104 


0.016 


0.090 




0.80 - 0.85 


0.11 


0.062 


0.013 


0.069 


0.85 - 0.90 


2.93 


0.025 


0.007 


0.039 


0.90 - 0.95 


0.76 


0.013 


0.005 


0.018 


0.95 - 1.00 


1.60 


0.012 


0.006 


0.005 



Table 5: As Table 1, but for the MARK II experiment [8] . The centre of mass energy is 29 
GeV and the overall normalisation factor from the fit is 1.044 (not included in the data 
column). The fourth column includes both the statistical and systematic errors. 
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18.93 
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152.6 
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1.06 


23.5 


0.4 


22.9 




0.20 - 0.30 


0.05 


7.79 


0.22 


7.77 


0.30 - 0.40 


0.02 


3.17 


0.15 


3.16 


0.40 - 0.50 


1.35 


1.23 


0.09 


1.34 


0.50 - 0.60 


2.10 


0.532 


0.050 


0.605 


0.60 - 0.70 


1.42 


0.310 


0.038 


0.267 


0.70 - 0.80 


1.40 


0.124 


0.022 


0.099 




0.80 - 0.90 


6.25 


0.021 


0.006 


0.036 


0.90 - 1.00 


5.65 


0.0038 


0.0018 


0.0081 



Table 6: As Table 1, but for the AMY experiment [16]. The centre of mass energy is 54 
GeV and the overall normalisation factor from the fit is 1.004 (not included in the data 
column). The fourth column includes both the statistical and systematic errors. 
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9.34 


262.0 


1.1 


2.8 


281.5 


0.03 - 0.04 


3.79 


184.3 


0.9 


1.4 


192.0 


0.04 - 0.05 


1.37 


136.7 


0.8 


0.9 


139.3 


0.05 - 0.06 


1.07 


103.0 


0.7 


0.6 


104.5 


0.06 - 0.07 


0.00 


83.3 


0.6 


r\ A 

0.4 


82.0 


0.07 - 0.08 


0.10 


67.7 


0.6 


0.6 


66.0 


0.08 - 0.09 


1.15 


56.1 


0.5 


0.3 


53.4 


0.09 - 0.10 


0.94 


47.0 


0.5 


0.2 


44.9 


0.10 - 0.12 


1.55 


37.0 


0.3 


0.2 


35.1 


0.12 - 0.14 


1.20 


27.9 


0.3 


0.1 


26.6 


0.14-0.16 


1.07 


21.3 


0.2 


0.1 


20.3 


0.16-0.18 


2.23 


16.8 


0.2 


0.1 


15.8 




0.18 - 0.20 


4.78 


13.71 


0.19 


0.13 


12.61 


0.20 - 0.25 


0.00 


8.93 


0.09 


0.12 


8.80 


0.25 - 0.30 


0.06 


5.43 


0.07 


0.08 


5.38 


0.30 - 0.40 


0.05 


2.88 


0.04 


0.04 


2.85 


0.40 - 0.50 


0.08 


1.24 


0.02 


0.02 


1.21 


0.50 - 0.60 


0.04 


0.534 


0.016 


0.012 


0.522 


0.60 - 0.70 


0.13 


0.230 


0.011 


0.004 


0.230 


0.70 - 0.80 


1.44 


0.090 


0.006 


0.002 


0.096 





Table 7: As Table 1, but for the ALEPH experiment [25]. The centre of mass energy is 
91.2 GeV and the overall normalisation factor from the fit is 0.984 (not included in the 
data column) . 



Source 


Error 


on 


Errors from Fit 


+20 
-11 


MeV 


Combinations of experiments 


±15 


MeV 


Heavy Quark Fractions 


±10 


MeV 


Fragmentation dependence 


±9 


MeV 


X dependence 


±10 


MeV 


Cut-off dependence 


±30 


MeV 


Scale dependence 


±40 


MeV 


Total 


±60 


MeV 



Table 8: Summary of systematic errors. The total error was obtained by adding quadrat- 
ically all errors. 
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Figure 1: a) The correction factor and b) the corrected inclusive momentum distribution 
where x = Phadron/ Ebeam from TASSO data at 35 GeV and DELPHI data at 91.2 
GeV. The sohd curves are results of the fits to the second order QCD matrix element. 
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DELPHI 




Figure 2: The ratio of the curves in Fig. lb. The dashed-dotted hue assumes that the 
flavour composition at 91.2 GeV is the same as the one at 35 GeV. As is apparent from the 
small difference between the solid and dashed-dotted lines, the increase in heavy quark 
production at the Z*^ resonance does not influence the scaling violation strongly. The 
reason is simply that although the heavy quark fragmentation is harder, the momentum 
spectra after the decays look similar to the ones from the light quarks and the difference 
does not show the characteristic energy dependence from the scaling violation. 
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Figure 3: The dependence of the inclusive momentum cross section in GeV^ for 
various x bins. For most data points the errors are smaller than the symbols. The solid 
curves are results of the fit to the data at high and high x, i.e. the data in the corner 
(bottom, right) indicated by the dashed lines. 
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Figure 4: The phase space for events with gluon radiation for two centre of mass energies 
(35 and 91 GeV) in the Xi versus X2 plane, where Xi and X2 are the fractional quark 
energies. The invariant mass between any pair of partons is required to be above 9.1 
GeV, which causes the empty bands near xi(2)—l- 
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Figure 5: Scaling violation is caused by the dependence of the gluon radiation. This 
figure shows the two contributions in first order: the increase in phase space of the qqG 
3-jet rate (solid line) and the running of the coupling constant, which causes a decrease 
of the 3-jet rate, if the fraction of phase space is kept constant (dashed line). The phase 
space was kept constant by requiring at all energies a minimum scaled invariant mass 
y = 0.08 between all partons. 



27 



M =MARKII 29 GeV A =AMY 54 GeV 

C =CELLO 35 GeV D ^DELPHI 91 GeV 

T -TASSO 35 GeV Al -ALEPH 91 GeV 



Exp, Comb. 


Result 


MTCADAI 


• 


MAD 


• 


CAD 


• 


MCAAI 


• 


MTCAD 


• 


CD 


• 


CTD 


• 


CTDAI 


• 


MTD 


• 



I I I I 

160 180 200 220 A In MeV 



Figure 6: The values obtained from fitting various combinations of experiments. 
The combinations were choosen such that the statistical errors were similar. They are 
indicated on the left using the abbreviations given at the top. 
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Figure 7: The renormalisation scale dependence of Ug (left) and the X^in dependence 
(right). Xmin is defined as the the minimum value of x used in the fit. 
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Figure 8: The phase space distributions of 4-jet events at an energy of 91 GeV for various 
values of the cut-off parameter y and the renormalisation scale factor /. 
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Figure 9: The ratio of the momentum spectra at energies of 500 and 91 GeV compared 
with the same ratio at energies of 91 and 35 GeV. 
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